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Abstract 


The stability of a single layer Newtonian fluid flow down an inclined plane lined with 
soft solid layer and pressure-driven two-layer plane Poiseuille flow throu^ a channel 
lined with soft solid layers are studied. The flow down an inclined plane system consists 
of a Newtonian fluid of density p , viscosity p and thickness R flowing adjacent to a 

soft solid layer of density p , modulus of elasticity E, viscosity and thickness HR 

which is fixed onto a inclined rigid plane having inclination angle of 6 . The linearized 
Navier-Stokes equations for the fluid and the elasticity equations for the soft solid layer 
are solved numerically, and the characteristic equation is obtained using the boundary 
conditions. The characteristic equation for the growth rate is a nonlinear equation, so 
analytical solutions cannot be obtained in general. Numerical solutions are obtained by a 
procedure using the solution from low wavenumber asymptotic analysis as the starting 
guess. It is found that soft solid layer has a stabilizing effect when the liquid flow down 
an inclined rigid plane exhibits instability. The stabilization is strongly affected by the 

nondimensional elasticity parameter T = , the nondimensional solid layer thickness 

ER 

p 

H, the ratio of viscosities r],=— and wavenumber k. 

The two-layer flow system consists of two layers of linear viscoelastic solid of 
thickness H^R and H^R , shear modulus E, and E ^ , and viscosity p^^ and p^^ fixed 

onto rigid surfaces at z = -iT,i? and z = (). + H^)R respectively, a layer of Newtonian 
fluid (fluid B) of thickness pR in the region Q<z<pR with viscosity p ^ , and another 
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Newtonian fluid (fluid A) of thickness (l-y0)i? in the region pR<z<R with 
viscosity . Similar to the previous analysis in this case also the Navier Stokes equations 

for the fluid layers and the elasticity equations for the solid layers are solved numerically, 
and the characteristic equation for each layer is obtained using the boundary conditions. 
The existence of interfacial instability in two-layer flow through rigid channel depends on 
viscosity ratio and thickness ratio (3 of fluid layers. Asymptotic results predict the 

soft solid layer has a stabilizing effect on the interfacial instability between two 
Newtonian fluids. Asymptotic results are valid only in the low wavenumber limit; we use 
a numerical procedure to continue the results to finite wavenumber. Numerical results 
reveal that two fluid interface can be completely stabilized by choosing appropriate 


values of parameter , solid layer thickness viscosities of solid layers 

E^R 


rf^ and interfacial tension between two-fluids E ((Z? = 1 and 2 are indices for solid layers). 
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Chapter 1 


Introduction 


Fluid flow adjacent to a flexible surface is often encountered in biological systems and 
biotechnological processes. The flow, of blood and other biological fluids in the body 
takes place through tubes whose walls are made of elastic materials such as tissues and 
membrane, and more recently, in microfluidic devices where soft elastomers are used in 
the fabrication of microchannels. A clear understanding of fluid flow past soft solid 
surfaces, and the way in which the solid interacts with the flow, will help in the accurate 
design and development of these applications. Such flows are also encoimtered in 
industrial applications which often involve fluid flow near polymer matrices or 
membranes. Recent experimental and theoretical studies indicate that the characteristics 
of the fluid flow in these systems could be very different from that in a rigid-walled 
channel due to the effect of the elasticity of the wall. In particular, the transition from 
laminar to turbulent flow in a flexible walled channel is qualitatively different from that 
in a rigid channel. A good understanding of the transition from laminar to turbulent flow 
in flexible channels could be important in technological applications because laminar and 
turbulent flows have very different transport coefficients. In applications where high heat 
or mass transfer rates are desired, it would be necessary to operate the system in the 
turbulent regime where the transfer rates are three to four orders of magnitude greater 
than that in the laminar regime. In processes where low drag forces are desirable the 
system could be operated in the laminar regime. The dynamics of fluid flow past soft 
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solids is very different from that of flow past rigid surfaces due to the deformability of 
the solid surface, which couples the dynamics of the fluid with the deformation of the 
solid. In particular, this coupling could lead to interfacial instabilities between the fluid 
and the solid, as first theoretically predicted by Ktunaran et a/. [5] and experimentally 
demonstrated by Kumaran and Muralikrishnan [6,9]. An interesting feature concerning 
these instabilities is that they can occur even in the creeping flow limit, in the absence of 
fluid and solid inertia. 

The problem of maintaining stable interfaces in two-fluid viscous flows is of 
considerable interest to many modem engineering flow processes. Interface stability is of 
critical importance in coextmsion processes in the plastics industry where multilayer 
flows are relevant. The coextrusion operation consists of combining several melt streams 
in a feedblock. The combined melt stream then flows to the die where the layers take 
their final dimensions. Under certain operating conditions, interfacial instabilities arise 
with detrimental effects on m.echanical, optical, and barrier properties of the final 
product. Efficient and robust control of the flow and suppression of instabilities is then of 
paramoimt importance in maintaining stable operating states. 

In the petroleum industry, for example, pumping costs in pipeline operations may 
be reduced considerably by adding a low viscosity fluid into the system. Experimental 
investigations of two-fluid pipe flow indicate that low viscosity fluids eventually 
encapsulate the more viscous fluids in both high and low Reynolds number flows [4], 
Interfacial instabilities in two-layer shear flows have been investigated by various 
authors. In his pioneering work, Yih [17] studied the stability of the interface using a 
long-wave asymptotic technique and showed that viscosity stratification alone can cause 
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instability even at vanishingly small Reynolds numbers. Subsequent to the study of Yih, a 
large number of studies (for an extensive review, see Joseph and Renardy [4]) have 
focused on the theoretical understanding of interfacial instabilities iii two-layer and 
multilayer flows of Newtonian and viscoelastic fluids. 

Free surface flow of fluids is of great importance in many industrial applications 
where solids are coated with polymeric films. For example, in slide coating operations 
the flow of liquid flow down an inclined plane is a fundamental step in production of 
multilayer products [2]. It is well-known that the rect i l i near flow of multiple fluid layers 
down an inclined plane can be unstable. The predominant application of the inclined 
plane flow geometry is in the manufacture of multilayer products such as photographic 
films, where the inclined plane provides a means by which individual fluid layers can be 
stacked on top of one another. At the end of the plane, these layers are then deposited 
simultaneously on a moving substrate [15]. 

In the present work, we consider gravity flow of a single layer Newtonian fluid 
down an inclined plane (shown in Fig 2.1) and pressure-driven two-layer plane Poiseuille 
Newtonian flow (shown in Fig 3.1). 

Linear stability of two superposed fluids in plane Poiseuille and Couette flow was 
first studied theoretically by Yih [17], who presented an analysis for two-dimensional, 
long-wavelength pertmrbations. Yih [17] derived a general expression for the complex 
wavespeed and showed that viscosity stratification can give rise to an interfacial mode 
depends on the geometric and physical parameters of the flow systems, notably the 
viscosity ratio and thickness ratio of fluid layers. In the analysis for two-layer 
Couette flow Yih [17] showed that even in tlie presence of a vanishingly small Reynolds 
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number, the interface becomes unstable in the long-wave limit due to viscosity 
stratification when the thickness of the more viscous fluid is smaller than the thickness of 
the less viscous fluid. The interface is stable when the thickness of the more viscous fluid 
is larger than that of the less viscous fluid. In the analysis for Poiseuille flow, Yih [17] 
did not discuss explicitly how the growth rate depends on the parameters p and jj .^ , and 
gave the numerical results only for the case fluids having same density and p = Q.5, 
thereby demonstrating effect of viscosity stratification alone is sufficient to cause 
instability. Yiantsios and Higgins [16] extended Yih’s [17] results for small wavenubers 
to large wavenumbers and accounted for differences in density and thickness ratios, as 
well as the effects of interfacial tension and gravity. Yiantsios and Higgins gave a 
relation between P and for a stable two-layer plane Poiseuille flow configuration. 

Papanastasiou et al. [10] performed the linear stability analysis of n-layer plane Poiseuille 
flow and had presented the asymptotic solutions for long and short wavelengths and 
numerical solutions for wavelengths of 0(1). Papanastasiou et al. [10] suggested that 
nonsymmetric arrangements, with thinner material in the outer layer and thicker material 
in the middle layer, tend to more unstable than symmetric arrangements of the same 
material. Recently Pinarbasi [10] presented the interface stabilization by surface heating 
and cooling. He did the analysis for two set of immiscible liquids silicone/water and 
oil/water (water at the bottom layer in both cases) and showed that an imposed wall 
temperature difference can be stabilizing or destabilizing depending on the disturbance 
wavenumber and layer thickness ratio. 
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The onset of instability of a liquid film flow down a stationary inclined plane was 
theoretically analyzed by Benjamin [1] and Yih [18]. In particular, they showed that the 
film on the vertical plate is always unstable and the Reynolds nmriber is greater than 

■^cot^ where 0 is inclination angle then ordy instability exists. When the inclined plane 
6 

has zero angle of inclination, the horizontal layer of liquid is stable. Lin et al. [7] carried 
out the stability analysis of viscous liquid film flow down an inclined plane that oscillates 
in the direction parallel to the flow. They showed that the onset of instability in a liquid 
film flowing down a stationary inclined plane can be suppressed by imparting an 
oscillation to the plane in a direction parallel to the flow. For a given liquid film, more 
than one band of frequency of oscillation may exist for the stabilization. Recently Lin et 
al. [8] carried out a similar analysis for two-layer film and demonstrated that the 
instability in a stratified film flow can be stabilized by imparting appropriate plate 
oscillation. 

Kumaran et al [5] studied the stability of the plane Couette flow of a Newtonian 
fluid past a deformable wall of finite thickness, which was modeled as an incompressible 
linear viscoelastic solid of fixed to a rigid substrate at zero Reynolds number. They 
showed that the plane Couette flow past a deformable solid is a finite wavenumber 
instability, and is stable for perturbations with very large and very small wavelengths at 
zero Reynolds number. Srivatsan and Kumaran [14] studied the same flow system but at 
nonzero Reynolds number. They used a numerical procedure and showed that the 
instability predicted by Kumaran et al. [5] at Re = 0 extends to finite and large Re. 
Shankar and Kumar [13] studied the stability of the two-layer plane Couette flow of 
Newtonian fluids past a deformable wall of finite thickness. They found that when the 
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two fluids are arranged so that they undergo low-A: instability in rigid channels due to 
viscosity stratification >1 >0.5 or fi^ <ly9&<0.5), the solid layer invariably 

has a stabilizing effect on two-fluid iuterfacial mode. When the two fluids are arranged 
such that the two-fluid interfacial mode is stable < 1 & > 0.5 or > 1 & < 0.5) , 

the solid layer could have both destabilizing and stabilizing effects, depending on its 
thickness. 

We want to explore the possibility of stability of stabilizing flow down an 
inclined plane, and two-layer pressure-driven Poiseuille flow using soft solid layers, 
following the previous work of Shankar and Kumar [13]. 

The outline of this thesis is as follows. In this work we first intend to study the 
problem of stability of liquid flow down an inclined plane lined with a soft solid layer in 
Chapter 2 and then in the next chapter we present the stability analysis of pressure-driven 
two-layer plane Poiseuille flow in a charmel lined with soft solid layers. Both the chapters 
consist of formulation of the linear stability analysis followed by presentation of 
asymptotic results and discussion of numerical method used to solve the governing 
stability equations. Results are then presented in successive sections primarily in tenns of 
plots between imaginary part of wavespeed and wavenumber and secondarily in terms of 
neutrally stability diagram. Conclusions are given in last section. 
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Chapter 2 

Stability of Liquid Flow down an Inclined Plane Lined 
with Soft Solid Layer 


In this chapter, we consider the stability of the Newtonian fluid flow past a soft, 
deformable solid layer which is fixed onto an inclined plane. Such a configuration might 
be relevant in the manufacturing of multilayer products such as photographic films [15]. 
The onset of instability of a liquid film flow down an inclined plane was analyzed by 
Benjamin [1] and Yih [18]. In particular, they showed that the film on the vertical plate is 
always unstable and when the inclined plane has zero angle of inclination, the horizontal 
layer of liquid is stable. In this chapter we are interested in studying how the instability of 
a liquid film flow down an inclined plane can be influenced by soft solid layer coatings. 
In particular, we explore the possibility of suppression of the instability by the soft solid 
layer. 

The organization of this chapter is as follows. The governing equations, base state 
and the linearized stability equations are presented in Sec. 2.1. In Sec. 2.2, we briefly 
discuss the numerical method used to solve the linearized stability equations, and in 
Sec. 2.3 we provide representative results for the complex wavespeed as a function of the 
wavenumber, and -secondarily in the form of neutral stability diagram. Conclusions are 
given in Sec. 2.4. 
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2.1. Problem Formulation 


The system configuration consists of a Newtonian fluid of density p , viscosity p 
and thickness R , flowing past a deformable wall (soft solid layer) having density p^ , 
viscosity , coefficient of elasticity E and thickness///?, fixed onto a inclined rigid 
plane at z = (l+H)R, in the region 0<z<R as shown in Figure 2.1. The liquid has a 
free surface at z = tj{x) . The soft solid layer is modeled as an incompressible linear 
viscoelastic solid, similar to that used in previous studies [5,13,14]. The soft solid wall is 
at rest in the unperturbed state but there is a strain in the soft solid due to fluid stress at 
the fluid-solid interface. The base velocity profile in the fluid is a parabolic function of 
z . Small perturbations to the flow are induced by spontaneous fluctuations, and we study 
the effect of wall deformability on the stability of these perturbations. 

2.1.1. Dimensional Governing Equations 

With the origin of Cartesian coordinates (shown in Fig. 2.1) at the free surface 
and with v,, pj. and characterizing the velocity field in the fluid, displacement 

field of viscoelastic solid, pressure field in fluid and the pressure in viscoelastic solid 
respectively, the dimensional governing mass and momentum conservation equations are 
as follows. 

Equation of continuity 

5,.V; = 0. (2.1) 

Equation of motion [Navier-Stokes equation] 
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p{d, +Vjdj)v. =-d,pj. +pd]v,+pg.. (2.2) 

The elastic medium (soft solid layer) is considered to be incompressible, so that the 
displacement field satisfies the following incompressibility condition 

0,«,=O. (2.3) 

Equation of motion [Linear viscoelastic model] 

Pg = -diPg + E d) u. + p^d,d^jU, + p^g, . (2.4) 

In this equation, the first term on the ri^t side is gradient of the pressure, the 

second term on the right side is the elastic stress due to strain, and the third term is the 
usual viscous stress due to gradient in the velocity. 

Boundary conditions used to solve these equations are as follows. At free surface 
z = r]{x) the normal and tangential stresses in the fluid must vanish i.e. 

tJ..nj = Q, (2.5) 

nXjnj-Q. ( 2 . 6 ) 

The viscoelastic solid is fixed to a rigid surface so, it satisfies the zero displacement 
condition at the rigid surface z = (l + iT) i? : 

M,=0. (2.7) 

The conditions at the interface (z = i?) are continuity of velocities and continuity of 
stresses. 


v.=5,«,, 

(2.8) 

(■ EjJ — tj n,y Hj , 

(2.9) 


(2.10) 
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where Ty and 11,^. are total stress tensor inside the fluid and viscoelastic solid respectively 
and these can be written as 


7!, = — n . <5.. H-r. , 

ij V y’ 

Vy=n (a,v. +a^.v,.), 

where ty is deviatoric stress tensor in fluid, and 

ay=-{E+fi^d,){d,Uj + dju,), 

where cXy is deviatoric stress tensor in soft solid layer. 


( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 


2.1.2. Base State 

With reference to Figure. 2.1, the base flow assumed steady, is parallel to x -axis, 
with a flat free surface the velocity varjdng only with z . Since the pressure gradient 
the X direction and velocity component parallel to z is zero, the Navier-Stokes equatio; 
are simply 

pgsm9-\-fx^^ = 0, ( 2.1 

dz 

^^-pgcosd. ( 2.1 

dz 

These equations can be integrated with the boundary conditions y^ = 0 at z = 0 a: 

— ^ = 0 at z = . The result is 

dz 

v^={pgs,in6llp){R^ -z^), (2.1 

and the average velocity is given by 
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v„ = pg P} sinO I'i/J,. 


(2.18) 


2.1.3. Nondimensional Governing Equations 

It is useful to nondimensionaiize the various physical quantities at the outset and 
the following scales are used for this purpose: R for length, average velocity for 

velocity, ^ for time and for stresses and pressure. This gives the following 

R 

nondimensional equations. From this point onwards, all the vjiriables are nondimensional 
unless it is explicitly stated that they are dimensional. 

Base velocity profile for fluid: 



v.=-(l-^'). 

(2.19) 


o 

II 

l>" 

(2.20) 

Equation of continuity for fluid and solid layer: 



II 

o 

(2.21) 


5,m, = 0. 

(2.22) 

Equation of motion: 



Fluid: 

Re(d, + v.dj )v, = -d^pj- + djV, + . 

(2.23). 

Soft solid layer: 


(2.24) 


The scaled stresses in fluid and soft solid are 

Ty =-PfS:j + (.d,Vj +djV,), (2.25) 
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(2.26) 


= -Pg + (p + Vr +5,»y) . 
where the various non-dimensional quantities are defined as: 


r = ^, Re = ^^ and T]r=— 
ER fj. fi 


Nondimensional boundary conditions are 
at z = 7 (x) [perturbed jfree surface] 


at z = \ + H [rigid surface] 


at z = 1 [fluid-solid interface] 


T =0, 

(2.27) 

r^=o, 

(2.28) 

u, =0, 

(2.29) 

u.=0. 

(2.30) 

i-i 

il 

(2.31) 

, 

(2.32) 

II 

(2.33) 

. 

(2.34) 


2.1.4. Linear Stability Analysis 

In this section we use a linear stability analysis to calculate the complex wave 
speed of perturbations to the interface between the elastic medium and fluid. In the 
analysis the velocity and displacement fields are expressed as a siun of mean and 
perturbation terms, i.e. 
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V,. =v.+v;, w, .=«,+«;. (2.35) 

All the perturbation quantities are expanded in the forms of Fourier modes in the 
X direction, and with an exponential dependence in time: 

v! = v.(z)exp[iA: (x— ct)], m! =«.(z)exp[zA: (x-ct)], (2.36) 

where k is wave number, c is wave speed and Vj(z)and m,.(z) are the eigenfunctions, 

determined by solving the conservation equations. Substituting Eqs. (2.35) and (2.36) in 
the conservation equation for the fluid, Eqs. (2.21) and (2.23), we obtain the linearzied 
equation for the fluid as: 

J.V.+ zA:v^=0, (2.37) 

Re{ik (v, -c)v^ + (i.vjvj = -ikPf + [d^ -k^]v^ , (2.38) 

Re{ik (v^ -c)k} = -dfi + [dl -F]k . (2.39) 

The above equations can be combined to give a single fourth-order, Orr-Sommerfeld-like 
equation forv_. : 

i k Re{{v^ -c) [d^ -k^]+ (d^y^) } V, = [d^ -ef K . (2.40) 

The linearized equations for the soft solid layer are obtained by inserting (2.35) and 
(2.36) in (2.22) and (2.24): 


dM.+ iku^=0. 

(2.41) 

-Re k^c%= 

(2.42) 

-Re ec\ = -dj^ + (^-77. ikc)[d^ -e]ii, . 

(2.43) 


These equations can be reduced to a single fourth-order differential equation for u . : 
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(2.44) 


Re k^c\dl -k^]u.^ +(" 7 , ikc)[dl -k^Yu.^ = 0 . 

The linearized equations for perturbations to stress fields in the fluid and soft solid layer 
are derived firom (2.25) and (2.26): 

A: V., f_ = -pf+2d.v., (2-45) 

= (^ -i kcTj^) (dM, +iku^), n_ = -p^ + (^ -i kcp^)! dM^ . (2.46) 

Now turning to the boundary conditions, the free surface condition must be 
applied at z = 7 not at z = 0 . These conditions are obtained by Taylor-expanding about 
the flat interface position in the base state. Because the gradient of shear stress of the 
pressiure in the base flow is not zero at z = 0 , additional terms appear as a result of the 
linearization about z = 0 . Hence (2.27) and (2.28) can be written in the forms 

(t U = -37 + «-V, +ikvY U= 0 , (2.47) 

(f_ )_, = -Pf L -,.0 -3 7 cot <9 + 2 d^v.^ -A:^ 2 7 = 0 , (2.48) 

where 7 the Fourier expansion coefficient for the free-surface position 
7 = 7 exp[iA:(x-ct)], and 2 = y/(//v„) is the nondimensional surface tension for fluid. 
The linearized kinematic condition at z = 0 (free surface) is given by 

ik[v^ (z = 0) - c ]7 = V. . (2 .49) 

Linearized boundary conditions in terms of perturbed quantities at z = (1 + 77) are 

M,=0, (2.50) 

w,=0. (2.51) 

Interfacial conditions for the perturbations at fluid-solid interface, z = h(x) Eqs. (2.31)- 
(2.34) reduces to 
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v.= —ikcu., (2.52) 

= -ikcu^, (2.53) 

t = n.> (2.54) 

71 = n^- (2.55) 


In these equations, the height of interface, h{x ) , is approximated by m. at z = l; 

this approximation is valid for small perturbations. The second term in left side of 
equation (2.53) represents the nontrivial contribution that arises as a result of the Taylor 
expansion of mean flow quantities about the unperturbed fluid-solid interface. This 
additional term is responsible for the instability of interface between the fluid and 
deformable solid wall [13]. 

The differential S 3 ^tem governing the stability problem consists of Eqs. (2.40), 
(2.44), (2.47), (2.48), (2.50), (2.51), (2.52), (2.53) and (2.55). It defines an eigenvalue 
problem for the complex wave speed c. The problem at hand is to determine c as a 
function of r,7j^,Re,'L,0,k and 77. It is not possible to solve this differential eigenvalue 
problem analytically. However, for long waves, it is possible to carryout an asymptotic 
analysis in the small parameter k , and this procedure yields analytical solution to c as a 
series in k The asymptotic analysis was carried out by Shankar [12] and this gives an 
anal)4ical expression for complex wave speed, but is valid only for long waves [k « l) ; 

c = c°+kc\ (2.56) 

where is purely real and c' is purely imaginary. From Shankar’s [12] analysis 

c“=3.0, (2.57) 

c'=/{(1.2i?e-cot^)-9r77}. (2.58) 
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The asymptotic results shows that in the ^ 0 limit, the parameters E and do not 

appear in c“ and c' . In the Eq. (2.58), first term proportional to Reynolds number and 
inclination angle 6 is inertial contribution while the second term proportional to TH is 
soft solid layer’s contribution. Stability of liquid flow down an inclined rigid plane 
depends upon the Reynolds number as well as on the inclination angle. The asymptotic 
results show that the soft solid layer’s contribution is always stabilizing in the limit of 
long waves. However, the long-wave analysis does not address the issue of whether the 
soft solid layer is stabilizing at all A: . To answer this question, we undertake a numerical 
study of the stability problem which is valid at all wavelengths. 

We use the asymptotic result as an initial guess for the numerical method through 
which we can continue up to finite values of wavenumber. The flow is unstable, neutrally 
stable, or stable according to whether c' is positive, zero or negative. 

2.2. Numerical Method 

We first briefly outline the method to solve the governing equations described in 
the previous section. There are two fourth-order differential equations for the fluid layer 
in the terms of v. and for solid layer in terms of u.. We use a fourth-order Runge-Kutta 

integrator with adaptive size control to obtain the solution. To solve a fourth-order 
differential equation we must know the function’s value and its first three derivatives at 
any particular value of independent variable. For the solid layer at z = l+Ff , we have 
from Eqs. (2.50) and (2.51) 

M^=0,J,M-=0- (2.59) 
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We use two different sets of higher derivatives dlv. and dlv. at z = \ + H , which 
will yield two linearly independent solutions to the displacement field in solid layer [14]: 

u,=0, dM, = 0 , d% = 1 , d% = 0 , (2.60) 

«, = 0 , rf.w. = 0 , d% = 0 , d% = 1 . (2.6 1) 

Using these conditions at z = 1+if , we use the Runge-Kutta method to integrate 
the ODE [Eq. (2.44)] in solid layer up to z = l. The displacement field in soft solid is 
then obtained as a linear combination of these two solutions. 


w. = A^uf^ + 4“^^ ’ 

(2.62) 

44 = 44.uf +44«f\ 

(2.63) 

44 = 444 *^ + A^dluf^ , 

(2.64) 

dlu. = A^dlifP + A^dluf^ . 

(2.65) 


Here 4 and are constants which have to be determined by the boundary conditions. 
Now we evaluate the velocity field of fluid v., and its higher derivatives, in terms of 
same constants 4 and 4 with the help ofboundary conditions at fluid-solid interface i.e. 
Eqs. (2.52)-(2.55). Using these two set of values at z = l, we integrate the 
ODE [Eq. (2.40)] in the fluid up to z = 0 and the velocity field v. in the fluid is then 

obtained as a linear combination of these two set of solutions. At z = 0 fluid velocity 
field must satisfy the boundary conditions (2.47) and (2.48). We can write the equations 
in matrix form as: 


^( 1 ) ^ ( 2 ) 
zx zx 

Y 0 ) Y 



( 2 . 66 ) 
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In this matrix, the subscripts 1 and 2 over the quantities refer to the solutions for 
the stress field in the liquid at gas-liquid interface calculated using the first and second set 
of boundary conditions respectively. For nontrivial solution of Aj andA^, determinant of 
matrix is set to be zero, which gives the characteristic equation. The characteristic 
equation is non-linear in c , so it is not possible to solve for c analytically. A Newton- 
Raphson iteration technique is used to solve the characteristic equation and obtain the 
complex wave speed, for specified values of F , Re, I,, 0, k and H . 

We use the low- A: asymptotic results as starting guess for the numerical 
procedure, and continue the low- k results numerically to finite values of k . We have 
compared the results from numerical code for the low values of k with Shankar’s [12] 
asymptotic results and foimd that numerical code is agreeing with the asymptotic results 
for small values of wavenumber. 

2.3. Results and Discussion 

It is instructive to briefly summarize the results of Yih [18] for the case of 
stability of liquid flow down in an inclined rigid plane. Yih’s [18] solution for long- 

waves shows that for the configuration where Re >— cot ^ long-wave disturbances will 

6 

be amplified. Also the free-surface flow down a vertical plane = 90°^ is unstable for 

all finite Reynolds numbers. Using the numerical procedure described above, it is 
possible to recover the case of flow down an inclined rigid channel, by setting the 
parameters T = 0 or // = 0 . Here, by setting T = 0 we achieve the conversion of the soft 
solid layer in to rigid solid layer, which has very high values of shear modulus (shear 
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modulus G oc — ) and i/ = 0 signifies the disappearance of soft solid layer, which is 


fixed onto the rigid surface at z = (l + . The results of the low wavenumber asymptotic 

analysis are used as starting guess for this numerical solution. In Figure 2.2 the imaginary 
part of wave speed is plotted as a function of wavenumber, for different values of 
Reynolds number. Figure 2.2 validates Yih’s [18] result and shows that the 


configurations is stable if Re <— cot 6 {Re = 0.1,0 = 45°) , while in opposite case where 

6 


Re >— cot 0 , exhibits instability for the low and intermediate k . Reynolds number has a 
6 


destabilizing effect at perturbations with long as well as intermediate wavelengths, as 

5 5 

shown in the same figure. In both cases (i?e <— cot^ and Re>—cot0) short waves are 

6 6 

stable even in the absence of gas-liquid interfacial tension. Figure 2.3 shows that free- 
surface flow down a vertical plate ( ^ = 90° ) is unstable for all finite Reynolds number. 
Similar to previous figure, here also increase in Re has a destabilizing effect at long and 
intermediate wavelengths and short waves are stable even with zero values of gas-liquid 
interfacial tension. Shankar and Kumar [13] showed that, in the case of two-layer flow 
only the presence of a sufficiently strong nonzero interfacial tension only can stabilize the 
short wave imstable modes. Variation of c,. with k for E 9^ 0 is presented in Figure 2.4. 

On comparison of Figure 2.4(a) (E = 0.l) with Figure 2.2 (E = 0), we find that the 

interfacial tension has no significant effect in the growth rate for the perturbations having 
long wavelengths, while in the case of short and intermediate waves it quantitatively 
affects the growth ratb. From Figure 2.4(b) it is more clear that interfacial tension 
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decreases the growth rate at high- A; but in this region the flow is already stable. So we 
can conclude that, for the present case interfacial tension does not qualitatively affect the 
instability of flow down an inclined plane. 

Now we examine the effect of soft solid layer on the interfacial instability df a 
liquid flow down an inclined plane. No work has been reported in literature related to free 
surface flow past a deformable wall. This effect was analyzed in the limit of long waves 
by Shankar [12]and the results are given in Eq. (2.56)-(2.58). That analysis showed that 
as k-^0, at leading order the wave speed is puerly real and is identical to Yih’s [1 8] 
result. The destabilizing effect of fluid inertia appears at 0{k ) . And the effect of the solid 
layer also appears at the same order, which is always stabilizing in the long-wave limit. 

In the present study we wish to examine whether the stabilization persists for 
intermediate wavelengths as well. Before we discuss the present results, it must be 
mentioned that shear flow past a soft solid layer in the creeping flow limit could become 
unstable when solid layer become sufficiently soft or when it is sufficiently tliick. 
Srivatsan and Kumaran [13] showed that this instability is not present for nondimensional 
thickness // < 1 . So, if we choose soft solid layers with nondimensional thickness H <\, 
we will not excite the instability of the interface between the fluid and the soft solid. 
Also, the instability of Kumaran et.al [5,14] exists only for T ~ 0(1) values. We will 
show below that for free surface flow the stabilizing effect of the solid layer occurs at 
value of r~O(10'-), which is much smaller than the critical value required to 

destabilizing the fluid-solid interface in case of shear flow. 

Here we present the cases where flow down an inclined rigid plane is xmstable and 
study the effect of the soft solid layer. We first discuss the effect of solid layer on the 
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liquid flow down a vertical plane and this is presented in Figure 2.5. Figure 2.5(a) and 
Figure 2.5(b) shows the variation of qwith k for ^ = 90° and i?e = 0.1 for different 

values of F . The soft solid layer is stabilizing only when the parameter F exceeds a 
critical value. In Figure 2.5(a), the solid layer with F = 0.001 has no significant effect on 
instability because here F is less than the critical value, while free surface is stable for 
F = 0.1 and 0.5 in low- A: limit. However value of c,. at AkI for F = 0.1 and 0.5 is 

larger than c. for the case when F = 0 . So the soft solid layer has destabilizing effect at 

finite- k , while short waves are always stable regardless of the parameter F . In Figure 
2.5(b), we show the effect of increasing F further. From this figure we conclude that 
complete stabilization (at all A ) of the free surface in the case of vertical plane can be 
achieved by making the solid layer sufficiently soft. Figure 2.6 shows that for 
stabilization it is not essential to have nondimensional thickness H >\ .Figure 2.6(a) and 
Figure 2.6(b) shows the variation of c.with A for // = 0.8, ^ = 45° and Re — X for 

different values of F . Figure 2.6(a) shows that while short waves are still stabilized, 
increase in F has a destabilizing effect on finite wave length fluctuations which is similar 
to the case of flow down a vertical plane discussed earlier. In the Figure 2.6(b), we show 
the effect of increasing F further. Indeed, compare with the rigid results show that the 
range of unstable A for F = 5 are stable in the case of rigid plane. Therefore, increase in 
F introduces new unstable modes for finite values of A . Therefore, if one were to choose 
the shear modulus of the solid layer, it is advantageous to have F<1 for the present 
configuration. On comparing Figure 2.5(b) and Figure 2.6(b) we find that in case of 
liquid flow down an inclined plane, high values of F has destabilizing effect at finite- A 
while there is no such destabilizing effect in case of vertical plane. Figure 2.7 shows the 
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variation of c,. with k for H = 2, 6 = 45® and Re = 1 for different values of F . In Figure 
2.6(a) for F = 0.5 c,. values at all k are negative while in Figure 2.7(b) c,. values are 

positive at finite- A:. So increase in H has a destabilizing effect on finite- A: modes. On 
comparing Figure 2.7(a) and Figure 2.7(b), we find that for the same value of = 1 , 
increase in F has a destabilizing effect on growth rate at finite- k - So, if we were to use 
soft solids for stabilizing the interfacial mode it is necessary to have nonzero F H , but 
there is an optimum value oiTH , and beyond this optimum value of Fi7 the soft solid 
layer has a destabilizing effect on the free surface at finite wavenumber as well as on the 
fluid-solid interfacial instability. Figure 2.8(a) and 2.8(b) shows that at Re = 5 the soft 
solid layer can not stabilize the instability at all wavenumbers; There is always there is a 
possibility of finite wavenumber instability. On comparing Figure 2.8(a) and Figure 
2.7(a) we can conclude that Re has destabilizing effect on finite-wavelengths 
fluctuations. This is further illustrated in Figure 2.9 where the parameter F = 0.1 is 
sufficient for stabilization of free surface at Re = l while it is less than the critical value 
at Re = 3 and 5. Nonzero value of gas-liquid interfacial tension has no significant effect 
on the instability at finite wavenumbers, as shown in Figure 2.10. 

It is useful to construct neutral stability curves demarcating stable and unstable 
regions in the T-k plane, for fixed values of ?j^,Re,'E,6mdH . Such plots allow us to 

select the parameter F (i.e. the shear modulus of the solid layer), for complete 
stabilization at all wavelengths. Neutral curves are provided in Figure 2. 11 -Figure 2.15. 
When F is increased beyond the neutral curve, there is transition from unstable to stable 
perturbations and an additional unstable region is identified when F is increased beyond 
a critical value for the finite wavenumber fluctuations. There are two neutral curves 
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separated by a stable region between them, where the flow is stabilized for perturbations 
with all wavelengths. The upper neutral curve corresponds to finite wavenumber 
instability that is absent in flow down a rigid plane. This finite wavenumber instability 
arises due to the soft sohd layer only. Gas-liquid interfacial tension does not qualitatively 
affect the neutral stabihty curve. However, there is a qualitative change in the neutral 
stability curve due to variation in the ratio of viscosities of the fluid and soft solid 
layer ( 77 ^). In Figure 2.11 where 77 ^ =0 the region where flow is completely stabilized 

ranges from T ~ 10“^ to F ~ 10"’ while in Figure 2.12 where 77 ^ = 0.5 it is from F - 10"^ 

to F-10°. So viscosity of the deformable wall is always stabilizing on the finite 
wavenumber instability of the free surface due to which the stable region between the two 
neutral curves increases with 77 ^. Figure 2.13 is for same set of parameters as in 

Figure 2.11 but for Re = 3 where the gap between the neutral curves has been reduced 
very much. This agrees with our previous conclusion that, Reynolds number has 
destabilizing effect on finite wave length fluctuations. Gas-liquid interfacial tension has 
no significant effect on neutral stability curve. This is illustrated in Figure 2.14. Figure 
2.15 is for same data set as in Figure 2.1 1 but for // = 3 , in this figure the upper neutral 
curve which is due to soft solid layer, starts appearing at F = 0.26 while for H = \ it 
appears after F = 0.44 . Hence increase in H also has destabilizing effect at finite wave 
length modes. It is also useful here to estimate the dimensional value of shear modulus 
that would be required to realize the present prediction. We consider the Figure 2.11 
where Re = 1 and 9 = 45® and the region where the modes are definitely stable ranges 
from F~10"^ to F~10"'. To estimate E we use p = l000kgm~^ and 
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fi-\Q kg mr^sec~' . By using i?e = 1 thickness of film comes out to be 0(l 0 ^ ) m and in 

the mentioned range of F, we conclude that E should be in the range 1 O'* -10^ Pa in 
order for all modes to be stable. Figure 2.16 and Figure 2.17 more clearly demonstrate the 
stabilizing effect of viscosity ratio rj^ and destabilizing effect of nondimensional 
thickness H at finite wavenumber instability. 

2 . 4 . Conclusion 

According to Yih’s [18] solution in the small wave number limit, liquid 
flow dovra an inclined rigid plane is unstable, even in the limit of small Reynolds 

number. The condition for instability is Re should be greater than— cot ^ , where 0 is the 

6 

angle with respect to horizontal. Our study concerning the stability of flow down an 
inclined plane lined with soft solid layer shows that it is possible stabilize this type of free 
surface flow instabilities by soft solid layer coatings. The asymptotic analysis for long 
waves shows that the deformable wall contribution is always stabilizing. The results of 
the low wave number asymptotic analysis were continued to finite and hi^er values of 
k using a numerical solution of the governing stability equations. We have found that, in 
the case of pertiurbations with finite wavelengths, the soft solid layer has destabilizing 
effect, while it is always stabilizing for the perturbation with very large wave lengths. 
The present configuration is always stable to short waves, regardless of soft solid layer 
coating. The nondimensional parameter F required to achieve the stabilization is 
typically very small compared to unity. The interfacial instability at the fluid-solid 
interface occur at F - 0(1) [5,14] therefore it is possible to choose a “band” of F value 


24 



in which fluid-solid interfacial instability is not excited. Surface tension of liquid-gas 
interface has negligible effect on the instability. Viscosity of soft solid layer is always 
stabilizing. Stability depends on the parameters, viz., T (elasticity parameter of soft 
solid), H thickness of soft solid layer and (soft solid layer’s viscosity parameter). 

Neutral stability curves were presented which help in tuning the solid properties so as to 
achieve complete stabilization at all wave-numbers. 



Figure 2.1. Schematic diagram of liquid flow down an inclined plane coated by soft solid 
layer (All the quantities are dimensional in the figure) 
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(b) vs k 

Figure 2.5 Instability in liquid flow down vertical plane lined by soft solid layer: 
variation of c. with k for H = 2,Re = OA,0 = 90^ ,1]^ =0 and E = 0. 








(b) c vs A: 

Figure 2.6 Instability in liquid flow down an inclined plane lined by soft solid layer: 
variation of c, with k for H=0.8(<1), /?e = l, 0 = 45°, 77 ^ =0 and Z = 0. 
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(a) c, vs k 



k 


(b) c, vs k 

Figure 2.7 Effect of solid layer deformability: c,. vs k for H = 2, Re = 1,0 = 45°, 77 ^ = 0 
and S = 0. 
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Figure 2.12 Neutral stability curves: F vs it for if - 2, i?e - 1, - 0.5, S - 0 and 

9 = 45° . Stable regions are denoted by S, unstable regions by U. 
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Figure 2a3 Neutral stability curves: T vs k for H =^2,Re = 3, ?j^=0,'Z = 0 and 
& = 45^ . Stable regions are denoted by S, unstable regions by U. 


2 = 0.0 
2 = 0.1 



Figure 2.14Neutral stability curves: F vs ^ for H = 2,Re-l,% -0 and 0^45 
Stable regions are denoted by S, unstable regions by TJ. 







Figure 2.15 Neutral stability curves: F vs k for = 3, ice = 1, 77 ^ = 0, E = 0 and 
9 = 45° . Stable regions are denoted by S, unstable regions by U. 


Figure 2.16 Effect of viscosity of soft solid layer on F vs A for Re = 1,9 = 45°, E = 0 
and H = 3. 







and 2 = 0. 
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Chapter 3 

Stability of Pressure-driven Two-Layer Plane Poiseuille 
Flow in a Channel Lined with Soft Solid Layer 


In this chapter, we consider the stability of pressure-driven two-layer plane Poiseuille 
flow of Newtonian fluids in a channel lined with soft solid layer. Multilayer flows of 
immiscible liquids are involved in several engineering application such as in multilayer 
extrusion of plastic films, in multilayer coating and in the transportation of oil. At certain 
processing conditions, wavy interfaces are observed in these multilayer flows and this is 
known as “interfacial instability”. In the present study we consider plane Poiseuille flow 
of two layers of Newtonian fluid having different physical properties and study the effect 
of soft solid layer on the interfacial instability. 

This chapter is stmctured as follows: The governing equations, base state and the 
linearized stability analysis are presented in Sec. 3.1. In Sec. 3.2, we present the low 
wavenumber asymptotic results. In Sec. 3.3, we briefly discuss the numerical method 
used to study finite wavenumber perturbations, and in Sec. 3.4 we provide representative 
results for the complex wavespeed as a function of the wavenumber, and neutral stability 
diagrams in appropriate parameter space. Conclusions are given in Sec. 3.5. 
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3.1 Problem Formulation 


The system under consideration (shown in Figure 3.1) consists of two layers of linear 
viscoelastic solid of thickness H^R and Ffji? , shear modulus and E 2 , and viscosity 

and fixed onto rigid surfaces at z = -H^R and z = (l+H 2 )i? respectively, a 
layer of Newtonian fluid (fluid 5) of thickness PR in the region 0 < z < /?i? with 
viscosity , and another Newtonian fluid (fluid A) of thickness (1-P)R in the region 
PR <z <R with viscosity . We have assumed that the densities of the fluids are equal 
( p„ = Pf, = p), in order to exclusively focus on tihe instability of the two-fluid interface 

due to viscosity stratification. The soft solid layer is modeled as an incompressible linear 
viscoelastic solid, similar to that used in the previous chapter as well as in the previous 
studies [13]. The base flow consists of two parabolic profiles. Similar to the previous 
problem, small perturbations to the flow are induced by spontaneous fluctuations, and we 
study the effect of walls deformability on the stability of these perturbations. 

3.1.1 Dimensional Governing Equations 

The governing equations for fluid and solid layer are same as in 
Chapter 2. 

Equation of continuity 

5,v“ = 0. (3.1) 

Equation of motion [Navier-Stokes equation] 

pid,+vpjX=-diP^f+p°d]v^. (3.2) 
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Here oc — a and b represents the upper and lower fluid respectively. The soft solid layer 

is considered to be incompressible, so that the displacement field satisfies the following 
incompressibility condition 

5,< = 0. (3.3) 

Equation of motion [Linear viscoelastic model] 

= -S^pI +E^d]uf +jup,d].uf . (3.4) 

Here <p = l and 2 are indices for soft solid layers at top and bottom plates respectively. In 
theses equations v°,uf,p“ andpj characterizes the velocity field in the fluid layers, 

displacement field of solid layers, pressure field in fluid and the pressure in viscoelastic 
solid layers respectively. The soft solid layer (1) is fixed onto a rigid surface and it 
satisfies the no slip condition at z = -H^R 

mW=0 (3.5) 

The continuity of velocities and stresses at fluid-fluid and fluid-solid interfaces give rise 
to following equations 

at fluid-solid interface (z = 0) continuity of velocities and stresses 

( 3 . 6 ) 

n(%,= 3;‘n, , (3.7) 

at fluid-fluid interface (z = j3R) continuity of velocities and stresses 

vf=vr (3.8) 

7tnj = T;n, (3.9) 

at fluid-solid interface (z = i?) continxiity of velocities and stresses 
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The soft solid layer (2) satisfies no slip condition at z = (1+772 )i2 

Mf^=0. (3.12) 

The dimensional total stress tensor in the two fluids is given by 

T-=-pJS.+r^, (3.13) 

(S.vJ + aX), (3.14) 

where r“ is the deviatoric stress tensor in the fluid a and is pressure in the fluid a . 
The nondimensional total stress tensor 11, “ in solid layers is given by 

n;=-^;<s„+<. (3-15) 

+ a^), (s-is) 

where cr? is the deviatoric stress tensor in the soft solid layer xp . 


3.1.2 Nondimensional Governing Equations 

The scales used to nondimensionalize various physical quantities are as follows: R for 
lengths and displacement, v,, = — for the velocities where G is constant pressure 


gradient, — for time, and for stresses and pressure. 77, and7r2> therefore, are 

V. R 


nondimensional thickness of solid layers, while (l-^) aod ^ are, respectively, the 



nondimensional thickness of fluid A and B . From this point onwards, all the variables 
are nondimensional unless it is explicitly stated that they are dimensional. 

The Nondimensional governing equations in the two fluids are, respectively, the 
Navier-Stokes continuity and momentum equation: 


W =0 

(3.17) 

Reid, +v* dX=-d,p‘+fi’:d)v‘ 

(3.18) 

Here, " is the Reynolds number based on viscosity of 

fluid B, and 

U u Llf 

= — = 1 are the viscosity ratios.The nondimensional total stress tensor 

Ty in the two fluids is given by 


ii 

1 

+ 

(3.19) 

r^=M:id,v^ +djvn, 

(3.20) 

The nondimensional governing equations for soft solid layers are same as in 

Chapter 2. The displacement field satisfies the continuity. 

d,uf=0 (3.21) 

The momentum conservation equation is given by 

Re(d^uf) = djni 

The nondimensional total stress tensor O? in solid layers is given by 

(3.22) 


(3.23) 

< =(3- + ’j’b.y.ax +3i“J) 

(3.24) 


v 
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Nondimensional boundary conditions are: 
at z = -//, [bottom rigid plate] 


wf*=0, Mj')=0 (3.25) 

at z = 0 [Fluid j5-Solid (1) interface] 


at z = /? [fluid /1-Fluid B interface] 


at z = 1 [Solid (2)-Fluid A interface] 


II 

(3.26) 


(3.27) 

n*2 

(3.28) 

II 

^1 

(3.29) 

v; = vf 

(3.30) 


(3.31) 

rpa rpb 

(3.32) 

n 

II 

(3.33) 

v:=d,u^j^ 

(3.34) 

v;=8,»» 

(3.35) 

C 

H 

(3.36) 

r = 

(3.37) 


at z = \ + H 2 [top rigid plate] 

i7f)=0, = 0 (3.38) 
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3.1.3 Base State 


The steady velocity profiles in the two fluids are simply the Poiseuille flow velocity 
profiles: 

v‘=-z'+(!>z, (3.39) 

(3-40) 


where <!> = 


^U-1) 


3.1.4 Linear Stability Analysis 

In this section, the stability problem is formulated by the method of normal modes 
[3]. The velocity and other d)mamical quantities in the two fluids and the solid layers are 
decomposed into two parts, tiie base fields and infinitesimal perturbations: 

vr=T+vr', (3.41) 

All the pierturbation quantities are expanded in the forms of Fourier modes in the 
X - direction, and with an exponential dependence in time: 

v" = v“(z)exp[i^ (x-ct)], wf =Mf(z)exp[iJ: (x-ct)], (3.42) 

where k is wave number, c is wave speed and v“(z)and m“(z) are the eigenfunctions, 
determined by solving the conservation equation. For simplicity only two dimensional 
perturbations are considered. Substituting Eqs. (3.41) and (3.42) in the conservation 
equation for the fluid, Eqs. (3.17) and (3.18), we get the linearzied equation for the fluid 
as: 

rf.v“+i,tv“=0 (3.43) 
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Re{ik (vf -c)v; + (t/,v;)vf } = -ikP^ +M“[dl -k'^K (3-44) 

Re{ik (v; -c)v;} = -d,P; +^^[dl -e]v^ (3.45) 

The above equations can be combined to give a single fourth-order, Orr-Sommerfeld-like 
equation for v. : 


i k Re{{v^ -c) [d] -^^]+ {d]v^) } vf =^‘^[dl -k^f (3.46) 

The linearized equations for the soft solid layer are obtained by inserting (3.41) and 
(3.42) in (3.21) and (3.22): 

dfil+ikul=^ (3.47) 

-Re ec^u^ = -ikpl + ikc)[dl -k^]u^^ (3.48) 

^ (p 

-Re k^c^m = -d,pl + (—-77^ ikc){dl -k^]u^ (3.49) 

f /n 


These equations can be reduced to a single fourth-order differential equation forw" : 


Re k^c^[d^-k 


'KH- 


-p'^ikc)[dt-k^Yul=Q 


(3.50) 


The linearized equations for perturbations to stress fields in the fluid and soft solid layer 
are derived from (3.19) and (3.23): 


K =p:{d^r,+ikr^), t: = -p;+ 2 /iXv^ 


(3.51) 


n: =(p-ac>7.')(«+iMn. ^l=-Pl*^p-ikcri:)2dX (3.52) 


The linearized boundary conditions at the unperturbed interface position z = 0 between 
fluid B and the solid layer B are obtained by Taylor expanding about the flat interface 
position in the base state. These become 
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(3.53) 


v*= -ikcu^\ 

x;+Kv‘Ui;» = -itcc<.'', 


fh _ ] q ( 1 ) 


(3.54) 

(3.55) 

(3.56) 


In these equations, the height of interface, h^ (;c) , is approximated by at z = 0 ; this 

approximation is valid for small perturbations. The second term in left side of equation 
(3.54) represents the nontrivial contribution that arises as a result of the Taylor expansion 
of mean flow quantities about the unperturbed fluid-solid interface. This additional term 
is responsible for the instability of interface between the fluid and deformable solid wall. 

Similarly, the linearized boundary conditions at the unperturbed interface position 
at z = /? between the two fluids A and B are given by 


5.’ =?;. 

(3.57) 


(3.58) 

II 

(3.59) 


(3.60) 


where g is the Fourier expansion coefficient for the interface position 
g = gexp[j7(:(x-c/)], and E = is the nondimensional interfacial tension between 

fluids A and B. The linearized kinematic condition at z = /? between the two fluids is 
given by 




(3.61) 
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The linearized boundary conditions at the unperturbed interface position at z = 1 between 
the fluid A and solid layer A, siniilar to (3.53)-(3.56), are given by 


a 

II 

1 

' 

(3.62) 

= -ikcuf^. 

(3.63) 

7;:=n«, 

(3.64) 


(3.65) 

The boundary conditions at z = (l -t- ) are simply 



(3.66) 

Wlrile the boundary conditions at z = -i7, are 

= 0, = 0. 

(3.67) 


Differential equations (3.46) and (3.50) along with boundary conditions (3.53)-(3.67) 
completely specify the stability problem of interest in this study. 

The complex wavespeed c is a function of r^,T]^,Re,'L,^,k,H^ and ju ^ . For arbitrary 

and k, there are no closed form solutions to the fourth-order differential equation (3.49), 
and so a numerical method must be used to solve the stability problem in general. 
However, when we consider very long waves, i.e., k«\, an asymptotic analysis in the 
small parameter k is possible, similar to the analysis of Yih [17], which yields an 
analytical expression for the wavespeed as an asymptotic series in k. In the following 
section, we briefly discuss the low wavenumber asymptotic results carried out by 
Shankar [12]. These low-A: results are then used as starting guesses for a complete 
numerical treatment of the stability problem. 
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3.2 Asymptotic Results 

The as 5 aiiptotic analysis carried out by Shankar [12] gives an analytical 
expression for complex wave speed, but is valid only for long-waves {k « l) ; 

c = c“ + fe', (3.68) 

where leading order c° is purely real and first correction c* is purely imaginary. The 
asymptotic analysis shows that at in the ^ 0 limit, the parameter E (interfacial 

tension) and do not appear in c° and c' . The analytical expression for c' is very 
complicated when remaining parameters T^,'n%Re,'L,p,k,H„and jj.^ are left 
unspecified. However, when only j3 and /j.^ are specified, it is possible to obtain 

anal 3 dical expressions for c‘ with the other parameters left unspecified. First we consider 
two layers of fluid with single soft solid layer fixed onto bottom plate. For 
and p = 0.6 , the results for c° and c’ are 

c® =0.168052, (3.69) 

c' =j{-0.00283 r,(iT, -1.2368)iT,(iT, +0.01 897) + 0.00001 792 Re ] . (3.70) 

The first correction c* , has two contributions; the term proportional to Re is destabilizing. 
However, there is another term proportional to F, and if, which appears due to the 
presence of the solid layer placed at bottom plate. This soft solid layer contribution gets 
change by changing the arrangement of solid layer. Since c° is independent from and 

and depends only on P and //,. , for same values of P and but solid layer at top 

plate value of c® is same. The analytical expression for first correction c is very different 
and is given by 
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c' -/|-0.01536r /^2-0.01119r2F2-0.00283ri/2^H-0.00001792i?e} (3.71) 

On comparing Eq. (3.70) with (3.71), we find that for p = 0.6 and //, = 2 if we place the 

solid layer on bottom plate then solid layer is stabilizing only when nondimensional 
thickness H is greater than 1.2368 while in the case of solid layer fixed onto top plate, it 
is always stabilizing irrespective of values of if. 

Now we consider the configuration where fj,^ = 0.25 (< l) . First we consider two 

layers of fluid with single soft solid layer fixed onto bottom plate. For 
= 0.25 and p = 0.4 , the results are 

c° =0.458211 (3.72) 

c‘ =i {-0.1701 r, if, -0.1574rF/-0.046ir/ff + 0.0008859 i?e} (3.73) 

For same values of P and but solid layer at top plate, the analytical expression for first 

correction c' is 

c' =i {-0.04610 T^{H^ +0.2176) +0.0008859 Re] (3.74) 

In this case where fj.^ is less than 1 it is advantageous to fix the solid layer onto bottom 

plate because it is stabilizing for all values off/,, while it is destabilizing if < 0.6318 

when it is placed witib the top plate. 

In general , therefore, the asymptotic results can be expressed as follows: 

c = ^+ki\T^H^F[H^)+ReC,], (3.75) 

where c° is purely real, F(ifp) is a real valued function of (which can be positive 
or negative depending upon ff^) that represents the effect of solid layer, and C, is a real 
constant that represents the contribution due to viscosity stratification between toe two 
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fluids. A few representative results for and C, are summarized for different 

values of (>1 as well as <1) and p in Tables 3.1 - 3.6. From these results (and from 
the results for other values of and p which are not presented here), the general trend 

appears to be that when the two fluid layers are arranged so that they undergo the long- 
wave two-layer viscosity stratification instability, the solid layer usually has a stabilizing 
effect in the long-wave limit. The stabilization due to presence of solid layer depends on 
placement of solid layer, if it is placed next to more viscous fluid then its thickness has no 
qualitative effect on stabilization, while if it is placed near to less viscous fluid then 
stabilization occurs only when nondimensional thickness exceeds a critical value and 

this critical vale dep)ends on thickness ratio if fluids {^0) . On the other hand, when the 

two fluid layers are arranged so that the interface is stable in the long-wave limit, then the 
solid layer usually has a destabilizing effect, but this destabilization is also affected by 
placement of solid layer. If solid layer is next to more viscous fluid then it is destabilizing 
irrespective of thickness in the long-wave limit, while it could be stabilizing for 

certain values of if it is placed near to less viscous fluid. 

For example, the result for =2 and y3 = 0.6(see Table 3.1) shows that the 
solid layer placed at top plate has a stabilizing effect for all values of while solid 
layer at bottom plate is stabilizing for /f, >1.2368. In the same table for 0 = 0.2, solid 
layer at bottom plate is destabilizing for all values of H^ and solid layer at top plate is 
stabilizing for 0.9922. 
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The asymptotic results when solid layers are placed on both the plates for 
= 2 and /? = 0.6 is obtained simply by adding the contributions from top and bottom 
solid layers. This is possible owing to the linearity of the problem. 

c" =0.168052 (3.76) 

-0.00283 r, (77. -1.2368)^r, (if, +0.01897) 

c' =T -0.01536r2 ^z-O-OlllPT^ if^ -0.00283 T^ Hi > (3.77) 

+0.00001792 Re 

In Eq. (3.77) first two negative terms corresponds to lower and upper soft solid layer 
contribution and third term is inertial contribution. This asymptotic result is obtained by 
simply the addition of all the contribution known in Eq. (3.70) and (3.71). 

The above asymptotic results, however, are applicable only in the low-A: limit. In 
order to determine whether the predicted stabili 2 ation extends to perturbations with finite 
and high-it, it is necessary to solve the governing stability equations and boundary 
conditions numerically. In addition, at finite-A; the interfacial mode at the 
fluid 5-soIid (1) and fluid .4-solid (2) interface can become unstable when elasticity 
parameter 1'^^ is increased beyond a critical value. In the following section, we use the 

numerical solution to determine regions in the T^ -k plane where the two-fluid interfacial 
modes is stable. 

3.3 Numerical Method 

The asymptotic results fails to predict the results when wavelength of a 
disturbance is of 0(l) . To obtain the stability results over a wide range of wavelengflis, a 
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numerical solution of Orr-Sonunerfeld equations is undertaken. We use the same method 
which has been used in the previous studies [13,14] as well as to solve the stability 
problem of liquid flow down an inclined plane. We use a fourth-order Runge— Kutta 
integrator with adaptive step size control to obtain numerical representations of the 
linearly independent solutions to the fourth-order ODEs. We recast the fourth order ODE 
in each layer as a system of first-order differential equations for the variables 

For solid layer (1), this system of equation can be integrated 

fromz = -i/, to z = 0 , provided the values of w^'^and its first three derivatives are 
known at z = . The soft solid layer B satisfies the boundary conditions 

= 0, and = 0 at z = -H ^ . We use two different sets of higher derivatives 

k 


and ^ at z = -H^ that are consistent with the boundary conditions at z = -if, , 

which will yield the two linearly independent solutions to the displacement field in soft 
solid layer B: 


= 0 , =0, = 1 , dlu^!'> = 0 , 

Z ' 44 . '44. '44 

(3.78) 

uP = 0 , d^P =0, d^dP = 0 , d^dP = 1 . 

(3.79) 

Using these conditions at z = -if, , we use the Runge-Kutta method to integrate the ODE 

[Eq. (3.50)] in solid layer R up to z = 0. The displacement field in 

obtained as a linear combination of these two solutions. 

soft solid is then 

dP=AdP'+A,dP“, 

(3.80) 

d^dP=A,djiP'+Mr«P‘'’ 

(3.81) 

d^dP^^AdjdP'+^d'uP", 

(3.82) 






(3.83) 


here A, and are constants which have to be determined by the boundary conditions. 
Now we evaluate the velocity field v* of fluid B, and its higher derivatives, in terms of 
same constants Aj and yf^with the help of boundary conditions at fluid 5-solid B 
interface i.e. Eqs. (3.53)-(3.56). Using these two set of values at z = 0, we integrate the 
ODE [Eq.(3.46)] in fluid 5 up to z = /? and the velocity field v* in the fluid is then 
obtained as a linear combination of these two set of solutions. 

The soft solid layer (2) displacement field satisfies 

= 0, and =0 at z = 1 + 7 / 2 . Therefore, we choose the following initial 

K 

conditions at z = 1 + 7/2 consistent with the boundary conditions for the displacement 
field 


=0, = 0, =1, dlii?'^ =0, 

(3.84) 

= 0 , = 0 , = 0 , = 1 . 

(3.85) 


‘ Using these conditions at z = 1 + 7 / 2 , we integrate the ODE [Eq. (3.50)] in solid layer A 
up to z = 1 . The displacement field in soft solid ^4 is then obtained as a linear combination 
of these two solutions. 

«f> = 4«f+4M?" ' (3.86) 

here and are constants which have to be determined by the boundary conditions. 
Now similar to previous steps we evaluate the velocity field of fluid A, and its higher 
derivatives, in terms of same constants A^ and with the help of boundary conditions 
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at fluid /4-solid A interface i.e. Eqs. (3.62)-(3.65). Using these two set of values at z = 1 , 
we integrate the ODE [Eq.(3.46)] in fluid /4 up to z = and the velocity field vf in the 
fluid is then obtained as a linear combination of these two set of solutions. 

At z = (5 fluids velocity field v“ and v* must satisfy the bormdary conditions 
(3.57)-(3.60). We can write the equations in matrix form as: 



-all 

V. 

V. 

vf 

vf +<.vf' 


nyall 

xz 

■^XZ 


fall 




(3.87) 


In this matrix, the subscripts I and II over the quantities refer to the solutions for the 
velocity and stress field at two-fluid interface calculated using the first and second set of 
boundary conditions respectively. For nontrivial solutions of and/44, 

determinant of matrix is set to be zero, which gives the characteristic equation. The 
characteristic equation is non-linear in c , so it is not possible to solve for c analytically. 
A Newton-Raphson iteration technique is used to solve the characteristic equation and 

obtain the complex wave speed, for specified values of r^r^, , Re, 2, >9, k, 


H, and 


3.4 Results and Discussion 

We first briefly recapitulate the results of Yih [17]. Linear stability of two 
superposed fluids in shear flow was first studied theoretically by Yih [17], who derived 
the general expression for the complex wavespeed in the low wavenxunber limit for two- 
layer Couette flow and showed the interfacial instability depends on geometrical and 
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physical parameters of flow system {0mdju^)md it exists if the thickness of more 
viscous liquid is smaller than that of the less viscous liquid 

(//^ > 1 & /? > 0.5 or //, < 1^ & < 0.5) . For the opposite case, viz., when the thickness of 
less viscous liquid is smaller than that of the more viscous liquid 


{/^r &. 0 > Q.5 ox >1&^<0.5), the interface is stable (“thin layer effect” [4]). 

These results are presented in Figure 3.2 in the form of neutral stability diagram. For 
plane Poiseuille flow, Yih [17] did not discuss explicitly how the interfacial instability 
depends on these parameters and gave numerical results only for the case fi = 0.S. 
Yiantsios and Higgins [16] extended the Yih’s [17] result for two-layer plane Poiseuille 

flow and showed that the interfacial mode is neutrally stable when u . Neutral 

0 ^ 


stability diagram (//, vs/?) for two-layer plane Poiseuille flow in rigid-walled channel is 
given in Figure 3.3. 

The results of the asymptotic analysis carried out by Shankar [12] are valid only 
in the low wavenumber limit. Using nxxmerical procedure described in the previous 
section, we obtain the results at finite wavenumber. The results of the low wavenumber 
asymptotic analysis are used as starting guess for these numerical solutions. First we 
represent the numerical results for flow in rigid-walled channel. Figure 3.4(a) and (b) 
shows that the interfacial instability exists in the case of flow through rigid-walled 
channel if the geometrical and physical parameters [P and//J lie in the unstable region 
shown in Figure 3.3. In Figure 3.4(a) and (b) interfacial instability at finite k is stabilized 
by nonzero E . In the following discussion we consider the parameters and //^) where 
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the interfacial instability exists in the flow through rigid-walled channel and we study the 
effect of soft solid layer on the interfacial instability. 

Now we present the results for the two-layer flow system with soft solid layer 
which is placed either on bottom plate or on top plate or on both the plates. In the Sec. 3.3 
we have aheady discussed how the placement of solid layer affects the instability. Effect 
of thickness of solid layer when it is placed only on the bottom plate is presented in 
Figure 3.5. For p = 0.4,//^ = 0.25 soft solid layer is stabilizing for all values of and is 

shown in Figure 3.5(a). For >9 = 0.6,//^ =2 soft solid layer is destabilizing in low 
wavenumber limit fori/, = 1.1 which agrees with the asymptotic results [Eq.(3.70)] while 
for i/, = 2 soft solid layer’s stabilizing contribution is enough to overcome the 

destabilizing inertial contribution. Now we shift tire position of solid layer from bottom 
plate to top plate and results are given in Figure 3.6. Here also we found agreement of 
numerical results with asymptotic results. In Figure 3.6(a) for >0 = 0.4,//^ =0.25 solid 

layer is destabilizing for = 0.5 while it stabilizes the interface in low wavenumber 
limit when i/j become 2.5. For /? = 0.6,//^ = 2 solid layer is stabilizing for all values of 
in the low wavenmnber limit. If we compare the results for the flow system with 

solid layer and without solid layer (rigid-walled channel) we find that solid layer 
stabilizes the interface in low wavenumber limit while it destabilizes the flow at finite-A. 

Now we move on to explore the effect of elasticity parameter on interfacial 

instability. Figure 3.7 and Figure 3.8 represents the results when solid layer is placed only 
on bottom plate and only on top plate respectively. Here we have considered sufficient 
thickness of solid layer so in all the cases c. values are negative in low wavenumber 
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limit. From both the figures we conclude that increase in parameter decreases the c,. 

in low wavenumber limit but it amplifies the distinbances at finite-lc. 

Now we discuss the flow systems with solid layer placed on both the plates 
Figure 3.9 simply present stabilized two-fluid interface in the low waveniunber limit by 
placing solid layers on the walls of channel where the interfacial instability exists in rigid 
channel. In Figure 3.10 we present the effect of interfacial tension on the instability when 
both of walls are deformable. Similar to previous studies [13] in the case of two-layer 
plane Poiseuille flow presence of strong nonzero interfacial tension can only stabilize the 
short wave instability. Effect of Reynolds number is given in Figure 3.11. Increase in 
Reynolds Number increases the growth rate at all k, thus it is always destabilizing. 
Figure 3.12 shows the effect of viscosity of solid layers on the instability for two different 
set of parameter. Qualitatively the viscosity of solid layers has no significant effect and 
quantitatively it has very little effect. So we can conclude that the viscosity of solid layers 
has no significant effect in interfacial instability. Up to now we have discussed the 
stability of two-fluid interfacial mode only there arte two more fluid-solid interfacial 
modes are present. Kumaran et al [5,14] showed that to excite theses modes elasticity 
parameter F should be 0(l) . In the present study we have considered very small values 

of parameter r , so here the other two modes will not be excited. 

It is useful to constract neutral stability curves demarcating stable and unstable 
regions in the T^-k plane, for fixed values of rj^ ,Re,'L, P snAH Such plots allow us 

to select the parameter F^^i.e. the shear modulus of the solid layer), for complete 
stabilization at all wavelengths. We present the neutral stability curve for the flow system 
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where both of the walls of channel are deformable. So we have two elasticity parameter 
r, and Fj and we are fixing the elasticity parameter of one solid to get F^ vs k plots. 
Neutral curves are provided in Figure S.lS-Figvure 3.15. The general trend appears to be 
when Fp is increased beyond the neutral curve, there is transition from unstable to stable 

perturbations. In Figure 3.13(a) for = 0.4, =0.25 when F, exceeds the value 
0.00067 the fluid-fluid interface becomes stable in the low wavenmnber limit from an 
unstable interface. Similarly in Figure 3.13(b) for P = =2 when Fj is increased 

beyond 0.00085 there is transition from unstable to stable perturbations but in the low 
wavenumber limit only. In both the figures perturbations at finite-A: are still unstable 
which can be stabilized only by the presence of interfacial tension and is presented in 
Figure 3.14. Effect of Reynolds number on neutral stability curve is shown in 
Figure 3.15. This figure shows that increase in Reynolds number increases the unstable 
region in F^ - A: plane. 

3.5 Conclusion 

In conclusion, the present study concerning the stability of two-layer plane 
Poiseuille flow through a channel where walls of the channel are lined by soft solid layer 
shows that solid layer has a profound effect on the interfacial mode due to viscosity 
stratification between the two fluids. Shankar’s [12] low wavenumber asymptotic 

analysis shows that the effect of solid layer appears at same order as the 

destabilizing effect due to fluid inertia in the asymptotic expression for the complex 
wavespeed c. Both these contributions could be stabilizing or destabilizing, depending on 
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the fluid viscosity ratio , relative thickness of the fluid layers P, and the parameters 
characterizing the solid layers . When the two fluids are arranged so that they 

imdergo a low-fc instability in rigid channels due to viscosity stratification, the solid layer 
invariably has a stabilizing effect on fluid-fluid interfacial mode, and in the present study 
we are interested only in the stabilizing behavior of solid layer. When the two fluids are 
arranged such that the two-fluid interfacial mode is stable, the solid layer could have both 
destabilizing and stabilizing effects, depending on its thickness. In the present study we 
explore only the stabilizing characteristics of solid layer. The placement of solid layer 
also affects the stabilization. When solid layer is placed next to more viscous fluid it 
exhibits stabilizing behavior irrespective of its thickness, while the placement of solid 
layer near to less viscous fluid can be destabilizing if its thickness is less than a critical 
value and this critical value depends on thickness ratio of fluids. The nondimensional 
parameter T^ required to achieve the stabilization of two-fluid interfacial mode is 
typically very small compared to unity in the limit of low k. Kumaran et <3/.[5,14] showed 
that when T^ ~ 0(l), the interfacial mode between fluid and solid layer also become 

unstable. In the present study we have considered T^ ~ so the interfacial mode 

between fluid and solid layer will not be excited. The results of the low wavenumber 
asymptotic analysis were continued to finite and hi^er values of k using a numerical 
solution of the governing stability equations. Numerical solutions showed that the solid 
layer has a destabilizing effect on two-fluid interfacial mode at finite and large k. These 
finite and short wavelength unstable perturbations can be stabilized only by the presence 
of a nonzero interfacial tension between the two fluids. Neutral stability curves were 
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TABLE 3. 1 Summary of results from low wavenumber analysis when solid layer is 
placed at bottom plate for = 2.0 


p 


c. 

0.1 

9.52x10-' +1.05x10-' /7,+ 3.59x10-' 

-1.66424x10^ 

0.2 

1.05x10-' +9.61x10-' i7,+ 2.61x10-' if,' 

-6.0140x10"' 

0.3 

5.39x10-' +3.33x10' if, + 5.93x10-^ if,' 

-8.0559x10"' 

0.4 

4.56x10-^ +5.02x10-' ii,+ 1.41x10-' if,' 

-1.6275x10"' 

0.5 

-3.25x10"' (-3.1328+ fi,) (-1.1172 +if,) 

1.1011x10"' 

0.6 

-2.84xl0-' (-1.2368+ if,) (0.01897 + if,) 

1.7927x10"' 

0.7 

-6.89x10-' (-0.6221 + if,) (0.3029 +ii,) 

1.3369x10"' 

0.8 

-7.79x10' (-0.3110 + /i,)(0.4153 +if,) 

6.6255x10"' 

0.9 

-3.52x10' (-0.1230 + ii,)(0.4692 +if,) 

2.9082x10“' 
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TABLE 3.2 Summary of results from low wavenumber analysis when solid layer is 
placed at top plate for = 2.0 


p 

F(H,) 

c. 

0.1 

3.59x10-' (-0.4079 + (0.4839 +77^) 

-1.66424x10^ 

0.2 

2.61x10-' (-0.9922 +/7j) (0.3208+772) 

-6.0140x10-" 

0.3 

5.94x10-^ (- 2.3457 + 772 ) (- 0.31912 + 772 ) 

-8.0559x10-" 

0.4 

1.4119 (-19.4078 + 772) (-13.1494 + 772) 

-1.6275x10-" 

0.5 

-3.24x10-^ ((3.625-1.3635 /)+ 772 )((3.625 +1.3635 

1.1011x10-" 

0.6 

-2.84x10-' ((2.10-0.98 i) + 772 )((2.1 0+0.98 

1.7927x10-" 

0.7 

-6.88x10-' ((1.65-0.83 i) + 772)((l.65+0.83 i) + 772) 

1.3369x10-" 

0.8 

-7.79 X 1 O ' ((1 .44 - 0.75 /) + 772 ) ((l .44 + 0-75 i) + ) 

6.6255x10^ 

0.9 

-3.53x1 0-' ((1 .33 - 0.70 r) + ) ((l .33 + 0.70 i) + ) 

2.9082x10^ 
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TABLE 3 3 Summary of results from low wavenumber analysis when solid layer is 
placed at bottom plate for = 0.25 


p 


c. 

0.1 

-0.01440 ((1.30-0.69/) +/7,)((l.30 + 0.69/) +77.) 

1.9614x10"" 

0.2 

-0.0463 ((1.37-0.73 /) +77,) ((1.37 + 0.73 /) +77,) 

4.0866x10"" 

0.3 

-0.0638 ((1.49-0.79/) + 77, ) ((1.49 + 0.79 /) +77,) 

5.5680x10^ 

0.4 

-0.0461 ((1.70-0.88/) +77,)((l.70 + 0.88 /) +77,) 

8.8859x10"" 

0.5 

-0.0151 ((2.19-1.05 /) +77,)((2.19 + 1.05/) +77,) 

9.2418x10"" 

0.6 

-9.90x10"' ((4.17 - 1.56 /) +77,)((4.17+ 1.56 /) +77,) 

3.6217x10"" 

0.7 

1.021x10"" (-7.8911+77,) (-3.7832 + 77,) 

-1.1989x10"" 

0.8 

4.55x10"' (-1.7711 +77,) (0.0878+77,) 

-1.9493x10"" 

0.9 

1.31x10"' (-0.6925+77,) (0.4732+77,) 

-6.8955x10"' 
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TABLE 3.4 Summary of results from low wavenumber analysis when solid layer is 
placed at top plate for = 0.25 


p 


c. 

0.1 

-1.44x10'' (-0.0643 + 7 / 2 ) (0.4614 +7/j) 

1.9614x10-" 

0.2 

-3.33x10'^ (-0.1673 + 7 / 2 ) ( 0 . 4157 - 1 - 772 ) 

4.0866x10-" 

0.3 

-6.38x10'' (- 0 . 3353 - 1 - 772 ) ( 0 . 3476 - 1 - 772 ) 

5.5680x10-" 

0.4 

-4.61x10'' (- 0 . 6319 - 1 - 772 ) ( 0.2176 + 772 ) 

8.8859x10-" 

0.5 

-1.52x10'' (-1.2696 + 772) (-0.1138 +772) 

9.2418x10"" 

0.6 

-9.90x10-" (- 3.6684 + 772 ) (- 1 . 6835 + 772 ) 

3.6217x10-" 

0.7 

5.72x10-' +1.49x10-' 772+ 1-02x10^ H\ 

-1.1989x10-" 

0.8 

2.83x10-' +2.13x10-' 772+4.54x10-' El 

-1.9493x10-" 

0.9 

4.08x10-' +4.22x10-' 772+ 1.31x10-' El 

-6.8955x10-' 
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TABLE 3.5 Summary of results from low wavenumber analysis when solid layer is 
placed at bottom plate for = 4.0 


p 



0.1 

1.02x10-^ +1.05x10-' 77,+ 3.28x10-' 77,' 

-1.07x10^ 

0.2 

7.07x10-' +5.33x10' 77,+ 1.14x10-' 77,' 

-3.04x10-^ 

0.3 

1.43x10-' +3.74x10"' 77,+ 2.55x1 0-'77,' 

-1.87x10'^ 

0.4 

-2.47x10^ (77, -3.6684)(77, -1.6835) 

5.87x10-" 

0.5 

-3.78x10-' (-1.2695+ 77,) (-0.1137 + 77,) 

1.44x10-' 

0.6 

-1.152x10-' (-0.6318+77,) (0.2176 + 77,) 

1.38x10-' 

0.7 

-1.59x10' (-0.3352 + 77,) (0.3476 + 77,) 

8.7x10-" 

0.8 

-1.15x10-' (-0.1674 + 77,)(0.4157 + 77,) 

6.38x10-" 

0.9 

-3.60x 10-' (-0.0643+ 77, )(0.4614 + 77, ) 

3.06x10-" 
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TABLE 3.6 Smnmary of results from low wavenumber analysis when solid layer is 
placed at top plate for fj.^ = 4.0 


p 



0.1 

3.28x10-" (-0.6925 + 7 / 2 ) (0.4731 +772) 

-1.07x10-" 

0.2 

1.11x10-" (- 1 . 7711 + 772 ) (0.0878 + 772) 

-3.04x10^ 

0.3 

2.55x10-" (- 7 . 8911 + 772 ) (-3.7832+772) 

-1.87x10^ 

0.4 

-3.78x10" ((4.1759-1.5624 z) + 772 ) ((4. 1759+1. 5624 1 ) + 772 ) 

5.87x10-" 

0.5 

-3.78x10-" ((2.1917-1.0524 z) + 772)((2.1917+ 1.0524 z) + 772) 

1.44x10-" 

0.6 

-1.11x10-' ((1.7071-0.8813 z)+772)((i.7071+0.8813 z)+772) 

1.38x10“" 

0.7 

-1.59x10-' ((1.4938-0.7919 :)+ 772 ) ((1.493 8 +0.79 19/) + 772 ) 

8.7x10-^ 

0.8 

-1 . 1 1 X 1 0-' ((1 .3758 - 0.73 54 z) + 772 ) ((1 .3758 + 0-7354 z) + ) 

6.38x10^ 

0.9 

-3.60x10-" ((1.3014-0.6945 z)+772)((l.3014+0.6945 z)+772) 

3.06x10-" 
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TABLE 3.7 ummary of results from low wavenumber analysis when solid layer is placed 
only at bottom plate for = 0.5 


p 


c. 

0.1 

-7.05x10 ' ((1.32 - 0.69 i) +H,)((130 + 0.69 z) +/f,) 

2.326x10“' 

0.2 

-1.55x10-' ((1.44-0.75 i) +i7,)((l.37 + 0.73 i) +H,) 

5.3x10“' 

0.3 

-1.37x10' ((1.65 -0.83 z) +iT,)((l.65 + 0.83 z) +77,) 

1.069x10“" 

0.4 

-5.678x10' ((2.10-0.98 z) + 77,) ((2.10 + 0.98 z) +77,) 

1.434x10“" 

0.5 

-6.49x10-" ((3.62 - 1.36 z) +77,)((3.62+ 1.36 z) +77,) 

8.88x10“' 

0.6 

2.82x10-^ (-19.4078+77,) (-13.1494 + 77,) 

-1.30x10“' 

0.7 

1.118x10“' (-2.3457 + 77,) (-0.3191 +77,) 

-6.44x10“' 

0.8 

5.23x10“' (-0.9922+77,) (0.3208+77,) 

-4.88x10“' 

0.9 

7.18x10“' (-0.4079+77,) (0.4839+77,) 

-1.33x10“' 
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TABLE 3.8 Summary of results from low wavenumber analysis 
placed only at top plate for = 0.5 


when solid layer is 
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Top Rigid Plate 


Fluid 


Fluid R 


Bottom Rigid Plate 


Top Rigid Plate 


Fluid .4 


Fluids 


Bottom Rigid Plate 


(b) Perturbed state 

Figure 3.1 Schematic diagram of two-layer plane Poiseuille flow in a channel lined with 
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Figure 3.2 Neutral stability diagram for two-layer Plane Couette flow for long- 
wavelength disturbances. Stable regions are denoted by S, unstable regions by U. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


p 

Figure 3.3 Neutral stability diagram for two-layer Plane Poiseuille flow for long- 
wavelength disturbances. Stable regions are denoted by S, unstable regions by U. 
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k 


(a) = 0.4, =0.25 and/?e=1.5. 



k 

(b) /3 = 0.6,fi^ =2 andi?e=1.5. 


Figure 3.4 Interfacial instability in the rigid channels: variation of imaginary part of 
wavespeed c with wavenumber k for two-layer Plane Poiseuille flow in rigid channels 

(r, =0andr, =0). 





Hi= 1.5 


(a) p = QA,u= 0.25, T, = 0.08, = 1 .5, S = 0, r 


Hi=0. 


Hj=2.0 




(b) = 0.6, /I, = 2, r, = 0.05, i?e = 1 .5, 1 = 0, r, = 0, if, = 0 . 

Figure 3.5 Effect of solid layer deformability on fluid-fluid interfacial mode: variation of 
c. with k for different values of thickness if, when the solid layer is placed only at 
bottom plate. 







(b) fi = 0.6,jii^=2,r2=0-OS,Re = l.5,I. = 0,r,=0,H,=0. 

Figure 3.6 Effect of solid layer deformability on fluid-fluid interfacial mode; variation of 
c, with k for different values of thickness when the solid layer is placed only at top 
plate. 







(b) ;0 = O.6,/^.=2,H,=2.O,i?e = l-5,2 = O,r,=O,//2=O- 

Figure 3.7 Effect of solid layer deformability on fluid-fluid interfacial mode: variation of 
c, with k for different values of parameter F, when the solid layer is placed only at 
bottom plate. - 
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(a) >0 = 0.4, = 0.25, = 2.5, i?e = 1 .5 ^ r, = 0, //, = 0 . 



Figure 3.8 Effect of solid layer deformability on flui' 


l-fluid interfacial mode: variation of 


c. 


with k for different values of thickness F, when ^ e so 1 ayer is p ace on y a op 


pi 








n =0.00, r2= 
ri=o.o 8 , r 


(a) P = 0.4, = 0.25, H,=2.5 = 0.5, = 1 .5, Z = 0, riP = 0, 77 ^^' = 0. . 


0.0002 


-0,0002 


^-- 0.0004 


-0.0006 


-0.0008 


- 0.001 


— r^ = 0X2=0. 

Fi = 0.05, = 0.08. 


(b) /? = 0.6, //, = 2, //, = 1 . 1, i /2 = 1 .5, i?e = 1 .5, 2 = 0, 77?’ = 0, 77^ ' = 0 . 

Figure 3.9 Effect of solid layer deformability on fluid-fluid interfacial mode: variation of 
c. with k when the solid layer is placed at both plates. 





(b) = 0.6, //, = 2, //, = 1 . 1 , = 1 .5, r, = 0.05, r, = 0.08, = 1 .5, t]):’ = 0, vr= 0 . 


Figure 3.10 Effect of interfacial tension when solid layers are placed on both the plates 








(a) = 0.4, /i, = 0.25,//, = 2.5, r, =0.08 = 0.5,r2=0.05,E = O,//^ = = 0 . 



Figure 3.1 1 Effect of Reynolds number when solid layers aie placed on bodi the plates. 
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(a) p = 0.4,//, = 0.25,//, = 2.5, r, =0.08 = 0,5,r2=0.05,i?e = 1.5, Z = 0 . 


0.0002 


- 0.0002 


^- 0.0004 


- 0.0006 


- 0.0008 


- 0.001 



(b) /? = 0.6,//, =2,//, =1.1 ,i /2 =1.5,r, =0.05, r, = 0 . 08 ,i?e = 1.5,Z = 0. 
Figure 3.12 Effect of viscosity of solid layers when solid layers are placed on both the 





(b) /? = 0.6, //, = 2, //, = 1 . 1 , = 1 .5, r, = 0.05, Re = 1 .5, E = 0, ri) ^ = 0, 7^ = 0 . 


Figure 3.13 Neutral stability curves for fluid-fluid interfacial mode: variation of F with 
k . Stable region is denoted by S, unstable region by U. 
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(a) = 0.4, yU, =0.25,//, =2.5, H^=0.5, T, =0.05,Re = 1.5,E = 0.05,?7i 


(h) J3 = 0.6, //, = 2, //, = 1 . 1 , = 1 .5, r, = 0.05, Re = 1 .5, Z = 0.05, 7]y=0,vy = ^- 


Figure 3.14 Effect of interfacial tension on neutral stability curves for fluid-fluid 
interfacial mode: variation of F with k . Stable regions is denoted by S, unstable region 
by U. 
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Chapter 4 


Conclusion and Scope for Future Work 


A linear stability analysis of single layer and two-layer flow of Newtonian fluid to two- 
dimensional perturbations has been conducted. The numerical solutions are computed for 
a wide range of wavenumber with the help of the initial guess from solutions of low 
wavenumber asymptotic analysis. 

The present study concerning the stability of gravity flow of a Newtonian fluid 
past a soft solid layer which is fixed onto a inclined plane shows that solid layer has a 
profound effect on the free surface flow instability. The soft solid layer has a stabilizing 
effect in the low wavenumber limit while it is destabilizing at finite-A: when the solid 
layer is sufficiently soft. Disturbances at large-A are stable regardless of soft solid layer 
coating. In the case of two-layer plane Poiseuille flow soft solid layer has a stabilizing 
effect when the rigid-walled channel exhibits the interfacial instability. 

In the present work we have considered the stability of Newtonian fluids. 
However in practical application, we may encounter the non-Newtonian fluids as well. 
For example, many biological fluids, fluids in industrial applications etc are viscoelastic. 
In this study we have considered single layer and two-layer of fluids but industries may 
deal with multilayer flow like in manufacturing of plastic films, manufacturmg of 
photographic films, multilayer extrusion, multilayer coatings etc. In all these applications 
interfacial instabilities are detrimental to product quality. Since soft solid layer has 
stabilizing effect on two-layer interfacial instability, it is instructive to perform the same 
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analysis for multilayer non-Newtonian fluids in order to extend the present results to 
more realistic system. 
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